Scaling behavior of the contact process in networks with long-range connections 
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We present simulation results for the contact process on regular, cubic networks that are composed 
of a one-dimensional lattice and a set of long edges with unbounded length. Networks with different 
sets of long edges are considered, that are characterized by different shortest-path dimensions and 
random-walk dimensions. We provide numerical evidence that an absorbing phase transition occurs 
at some finite value of the infection rate and the corresponding dynamical critical exponents depend 
on the underlying network. Furthermore, the time-dependent quantities exhibit log-periodic oscilla- 
tions in agreement with the discrete scale invariance of the networks. In case of spreading from an 
initial active seed, the critical exponents are found to depend on the location of the initial seed and 
break the hyper-scaling law of the directed percolation universality class due to the inhomogeneity 
of the networks. However, if the cluster spreading quantities are averaged over initial sites the 
hyper-scaling law is restored. 

PACS numbers: 



I. INTRODUCTION 



It is known both for equilibrium and nonequilibrium 
systems that the presence of long-range interactions leads 
to different critical behaviors compared to universality 
classes characteristic for systems with short-range inter- 
actions p|. This has been demonstrated for a paradig- 
matic model exhibiting a phase transition to a single ab- 
sorbing state, the contact process (CP), that has been 
introduced to model epidemic spreading without immu- 
nization 0, L*i this simple model, lattice sites have 
two states (active or inactive) and active sites may be- 
come inactive or may render neighboring inactive sites 
active. The absorbing phase and the active phase of the 
systems are separated by a continuous phase transition 
with the universal behavior of the directed percolation 
(DP) [J, la la]- Long-range interactions have been in- 
cluded by allowing the activation process for far-away 
inactive sites with probabilities that decay algebraically 
with the distance. The critical exponents of the corre- 
sponding absorbing phase transition have been found to 
depend continuously on the exponent controlling the de- 
cay of activation probabilities 0, [1] • 

An alternative way of realizing long-range interactions 
is when the dynamical process is defined on a network 
with long links that connect distant sites. Networks 
of this type are the scale-free networks constructed by 
preferential attachment pjj, where the critical behav- 
ior of CP is controlled by the degree-distribution [{J [l(| • 
Other models of networks are those which are composed 
of a (i-dimensional regular lattice and additional long 
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edges. These arise e.g. in sociophysics [TT| or in the 
context of conductive properties of linear polymers with 
crosslinks that connect remote monomers [121 ] . In gen- 
eral, a pair of nodes separated by the distance / is con- 
nected by an edge with a probabili ty p i ~ (3l~ s for large 

i dS> H^L HB> HH> E3, HM> IS> I21L [2^, 12^, 12^- we 

mention that the case s = corresponds to the Watts- 
Strogatz graph [25[, that displays the small- world phe- 
nomenon, although that model is constructed by rewiring 
edges rather than adding new ones therefore the result- 
ing graph may be disconnected. An intriguing property 
of these graphs for d = 1 is that in the marginal case 
s = 2, the intrinsic properties show power-law behavior 
and the corresponding exponents vary continuously with 
the prefactor (3. Indeed, this has been conjectured for the 
diameter flasa function of the number of nodes L, that 
means D(L) ~ L diain where the dimension d m - m depends 
on (3 [13]. Later, power-law bounds have been established 
for D(L) (2(|. For a class of cubic networks with s = 2, 
the algebraic grow th of the diameter has been explicitly 
demonstrated [24| . Moreover, the mean-square displace- 
ment of random walks in such networks has been found 
to grow algebraically in time with an anomalous random 
walk dimension that is characteristic for the underlying 
network. This behavior contrasts with Levy-flights in 
the respect that, here, the decay exponent s does not ex- 
clusively determine the diffusion exponent but the latter 
depends also on the details of the structure of networks if 
s = 2. As opposed to random walks, much less is known 
for interacting many-particle systems on such networks 
[23| . In particular, the behavior of nonequilibrium sys- 
tems possessing an absorbing phase transition has not 
cleared up yet. The aim of the present work is to investi- 
gate the contact process on these networks. On the basis 
of the scaling of diameter and mean-square displacement 
of random walks, we expect a nonequilibrium system pos- 
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sessing an absorbing phase transition, such as the contact 
process, to be characterized by altered critical exponents 
when defined on such networks compared to the corre- 
sponding one-dimensional model. We will demonstrate 
by Monte Carlo simulations that this is indeed the case 
for a class of regular networks which, concerning the di- 
ameter and random walks, are known to be described by 
power-laws, similar to s — 2 random networks. The ad- 
vantage of studying regular networks is that, at least, the 
diameter exponent and the random-walk dimension are 
exactly known here and unlike for random networks no 
disorder (sample) average is needed to carry out here. 

The rest of the paper is organized as follows. The 
networks to be investigated are defined in sectionQj] The 
model and the studied quantities are specified in section 
IIIII Results of numerical simulations are presented in 
section HVl and discussed in section IVl 

II. NETWORKS WITH LONG LINKS 

A. Definition of networks 

In this section, we specify the networks on which the 
contact process is studied. First, a regular, one dimen- 
sional lattice (periodic or open) with L sites is considered, 
where the lattice sites are numbered consecutively from 
1 to L. To this lattice, where the degree of sites, i.e. the 
number of edges emanating from a site, is 2, links are 
added one by one until all sites become of degree 3. Sites 
of degree 2 will be called in brief free sites, and k will 
denote a fixed positive integer. When constructing the 
networks, in general, pairs of free sites that have k — 1 free 
sites between them are connected iteratively [24[. In the 
following, the steps of this procedure will be described in 
detail. 

All the networks that we study are defined by means 
of aperiodic sequences, therefore, we start with a brief 
introduction of them. The aperiodic sequences that we 
need are generated by substitution rules on the letters of 
a finite alphabet A = {a, b, c, . . . } that assigns a word 
(a finite string of letters) w a to each letter a £ A. A 
(finite) aperiodic sequence is obtained by applying the 
inflation rule iteratively (finitely many times) starting 
with a single letter (which is letter b by convention). We 
use the two-letter sequence defined by the inflation rule 

J a — > w a = a6a(6a)™~ 1 , . 

° n : { b -> w b = a(ba) n ~\ { > 

where n is a positive integer and a three-letter sequence 
called tripling sequence generated by 

{a — > w a = aba 
b — > Wb = cbc (2) 
c — > w c = abc. 

The two-letter inflation rule with n = 1 generates 
the silver-mean sequence with the first few iterations 



b,a,aba,abaaaba,abaaabaabaabaaaba, etc, whereas with 
the choice n — 2, the well-known Fibonacci sequence is 
generated. 

Having defined the aperiodic sequences, we return to 
the construction of networks. First, the class of networks 
is defined where k = 1, that means, neighboring free 
sites are connected recursively. Here, a finite aperiodic 
sequence is chosen and the edges of the (finite) initial 
one-dimensional lattice are labeled consecutively by the 
letters of this sequence. Assume that all the edges which 
are labeled by letter a G A have a common length, de- 
noted by l a . Furthermore, let the edge lengths l a be 
ordered for the case of two-letter sequences as lb < l a 
while for the tripling sequence as lb < l c < l a . Now, let 
us find the closest pair of free sites (or, to be precise, 
one of the pairs with the shortest spacing between them) 
and connect them with an additional edge. This step is 
then iterated until all sites become of degree 3. It is clear 
that initially there is a multitude of pairs which are sep- 
arated by the shortest distance (lb) and, in the first few 
steps, these pairs are connected subsequently. Obviously, 
as the procedure goes on, the minimal spacing increases 
and longer and longer edges will form. 

We have studied the contact process on three different 
networks with k = 1: the silver- mean network, the Fi- 
bonacci network and the tripling network, which are con- 
structed by using the corresponding aperiodic sequence. 
The structure of these networks is illustrated in Fig. [T] 




abaaaba abaabaa aba 




a babaabaa baba 




a baebcabaa bccbcabca baebcaba 



FIG. 1: Finite silver-mean, Fibonacci and tripling networks 
(from top to bottom). 

In addition to this, we have considered networks with 
k = 2, as well. The k = 2 tripling network and the 
k = 2 silver-mean network are constructed as follows 
[24| . First, the sites of a one-dimensional lattice are la- 
beled with the letters of the corresponding sequence. The 
sites are grouped into blocks corresponding to words w a 
in the inflation rule. Then, sites belonging to one-letter 
blocks are renamed according to the reversed inflation 
rule w a — > a, where w a is the one-letter word corre- 
sponding to the block. In blocks composed of three sites, 
the two lateral sites are connected, and the middle one 
is renamed again according to the reversed inflation rule 
w a — > a, where w a is the word corresponding to the 
block. The above step is then iterated until only one 
free site (the central one) is left. The third network with 
k = 2 that will be investigated is the cubic Hanoi-tower 
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aba cbc aba abc cbc abc abacbc aba 




abaaaba abaabaa aba 




abababa babababa 



FIG. 2: Finite k — 2 tripling, silver-mean and Hanoi-tower 
networks (from top to bottom). 



network |23j that can be constructed following the above 
procedure with the inflation rule 

/ a -> w a = aba . . 

aH: { b->w b = b. W 

The three networks with k — 2 are illustrated in Fig. [21 



B. Diameter and random- walk dimension 

In the rest of this section, we shall survey some intrinsic 
properties of the above networks that are exactly known 
and are relevant with respect to the off-critical dynamical 
behavior of the contact process. 

Beside the distance x measured on the underlying one- 
dimensional lattice, another metric is the shortest-path 
length I between two sites which is the minimal number 
of links that have to be traversed when going from one 
site to the other one. The average length of the short- 
est path between two sites separated by the distance x 
scales in these networks as I(t) ~ x dmia where d m i n is the 
shortest-path dimension of the network. The diameter 
D(L) of a typical finite graph with L sites, which is the 
largest shortest-path length between any two sites grows 
also algebraically as D(L) ~ L'* mln . The average number 
of sites V (£) that can be reached in at most £ steps start- 
ing from a given site scales as V(£) ~ £ dg where d g is the 
graph dimension that is related to <i m i n as d g = 1 / c? m i n ■ 

The other property that we need is the random-walk 
dimension of the network. Let us consider a continu- 
ous time random walk on (infinite) networks, where the 
walker can jump with unit rate to any of the sites con- 
nected with the site it resides. The random-walk di- 
mension d w is defined through the asymptotical relation 
[(a; 2 (t))]t yp ~ t 2 / d ™, where x(t) denotes the displacement 
of the walker at time t and [(x 2 (t))] t yp = expln(a; 2 (t)) is 
the "typical value" of (x 2 (t)). Here, (•) denotes the aver- 
age over different stochastic histories for a fixed starting 
position, while the over-bar stands for the average over 
starting positions. Note that the expected value (x 2 (t)) 
does not exist if t > since the expected value of edge 
lengths is infinite (in infinite networks). This accounts for 



that the average of ln(x 2 (t)) is considered instead. No- 
tice that d TW is a dynamical exponent that relates time 
and length scale of random walks. 

The dimensions d m i n and d rw of the networks defined in 
the previous section have been exactly calculated [HI, [24| ; 
the corresponding numerical values are shown in Table |T] 
and Table HH 

III. CONTACT PROCESS ON NETWORKS 

A. Definition of the model 

The contact process is one of the earliest and simplest 
lattice models that belongs to the DP universality class. 
It is a continuous time Markov process on a state space 
{0, 1} S , where S is a finite or countable graph, usually 
Z d . A site with state 0(1) is called inactive(active) or, in 
the context of epidemics, healthy (infected). The dynam- 
ics is defined by nearest-neighbor transitions that occur 
independently with given rates. In d dimension, an in- 
fected site can spontaneously become healthy (1 — > 0) 
with rate 1 or can infect one of its neighbors (0 — > 1) 
with rate A/ (2d) (see Fig. [3]). This process is defined on 
the cubic networks under consideration in the way that 
infected sites are healed with rate 1 as before, whereas 
each nearest- neighbor site is infected with rate A/3, such 
that (for a site of degree 3) the total infection rate is A. 
In numerical simulations, this process is realized by ran- 

• o o o o • • o • • t 

• • O O • • ••• O t+dt 

m in x i 

FIG. 3: Allowed transitions and the corresponding rates in the 
one-dimensional contact process. Full(empty) circles symbol- 
ize active(inactive) sites. 

dom sequential updates: a randomly chosen infected site 
either becomes healthy with probability 1/(1 + A) or a 
neighboring site is attempted to be infected with proba- 
bility A/[3(l + A)]. For sites of degree 2 (the central site 
in k = 2 networks and the surface sites in open k = 2 
networks) and for sites of degree 1 (the surface sites in 
k = 1 networks), no update is attempted with probability 
A/[3(l + A)] and 2A/[3(1 + A)], respectively. 

Throughout the paper we measure the time in units of 
Monte Carlo (MC) steps that corresponds to N update 
attempts, where N is the actual number of particles at 
the beginning of the MC step. 

B. Studied quantities 

The quantities, the time-dependence of which we have 
measured at criticality, are the average number of ac- 
tive sites N(t) and the survival probability P(t), which 
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is the probability that there is at least one active site 
at time t. Furthermore, in the case of seed simula- 
tions, when initially a single site denoted by to is ac- 
tive and all other sites are inactive, we have also mea- 
sured the second moment of the distance of the grow- 
ing cluster with respect to the origin to, the so called 
spread R 2 n (t). To be precise, this quantity is defined as 
R 2 n {t) = (J2f=i ni{t)x1{t))/N{t), where x t is the distance 
between site i and the origin to, whereas is a binary 
variable which is one(zero) for active (inactive) sites. In 
the case of a single initial seed, the above quantities are 
expected to follow power-laws asymptotically: 

N m (t)~t^ (4) 
P m {t)~ t - 5 ™ (5) 
i&(t)~t a/ "». (6) 

In an inhomogeneous system the critical exponents rj mi 
S m and z m may be different for different initial active 
sites to. We shall see, that this is indeed the case, at 
least for the former two exponents. We have probed three 
different positions for the initial seed. First, it has been 
located at the central site(s), i.e. site (L + l)/2 for k = 2 
networks, where the number of sites L is odd, and sites 
L/2 and L/2 + 1 for k = 1, where L is even. We shall 
refer to this arrangement of the initial seed by the index 
'0'. Second, the initial seed has also been located at the 
site from which the longest edge of the (finite) network 
emanates. (Note that the edge connecting site 1 and site 
L in Fig. [T]is not the longest one since they are neighbors 
on the ring.) This initialization is indexed by T. Note 
that in this case, as the process starts, the spread jumps 
immediately to 0(x 2 ), where x is the length of the longest 
edge. To ensure the smooth increase of the spread, we 
have modified its definition so that Xi is the minimum 
of the distances measured from the two sites which are 
connected by the longest edge. Third, we have considered 
networks built on open one-dimensional chains, where 
the end sites are of degree 1 for k = 1 and of degree 2 
for k — 2. In these networks, the initial seed has been 
located at the surface, i.e. at site 1. The exponents 
corresponding to this arrangement are indexed by 's'. 

In addition to this, we have measured N(t) and P(t), 
the number of active sites and the survival probabil- 
ity, respectively, that are averaged over seed simulations 
started from all possible initial sites 1, 2, . . . , L. The cor- 
responding exponents are denoted by r] av and <5 av . Note 
that the average of R 2 {t) over all possible initial sites 
diverges in an infinite network for any t > since the ex- 
pected value of edge-lengths is infinite (24j ■ Nevertheless, 
the averaging would not provide any new information 
anyway on the spread since, according to the numerical 
results, the dynamical exponent z m = z is independent 
of the location of the initial seed. 

Another dynamical scaling exponent characterizes the 
critical system that is started from a homogeneous, fully 
occupied initial state. In this case, the density pit) of 



infected sites decays asymptotically as 

p(t) ~ t~<*. (7) 

In case of models of the DP class defined on regular lat- 
tices, a = S holds due to the rapidity reversal symmetry 
(see [l|). We shall, however, see that this equality does 
not hold in general for the networks under study. 



IV. RESULTS 
A. Off-critical behavior 

According to our numerical results, below the criti- 
cal value of the activation rate A c , that is characteristic 
for the underlying network, the system is in the inac- 
tive phase. Here, the number of active sites is found 
to decrease exponentially in time similar to regular lat- 
tices. The dynamics in this phase can be essentially de- 
scribed by random walks of the infection since active 
sites typically become rapidly inactive after activating 
a neighboring site. Accordingly, the spread is well ap- 
proximated by the mean-square displacement of random 
walks, R 2 (t) ~ £ 2 / d ™, with the only difference compared 
to regular lattices is that the dimension d ivr entering the 
above relation is the anomalous random- walk dimension 
of the underlying network. 

In the active phase, A > A c , the inactivation processes 
are irrelevant and the infection is spreading with a con- 
stant speed across the network. Since in t time steps all 
sites within the distance I ~ t are activated, the number 
of active sites grows in time as N(t) ~ V(t) ~ fV^in. 
As the growing cluster of active sites is compact in this 
phase, the spread increases in time as R 2 (t) ~ f 2 / d >m«n. 

The above laws for R 2 (t) below and above the critical 
point provide the bounds d m i n < z < d lw for the non- 
trivial critical dynamical exponent z. Furthermore, the 
number of active sites in surviving samples must not grow 
faster at criticality than in the active phase, which yields 
the inequality rj + S < l/d min . We shall see, that the 
measured critical exponents are compatible with these 
(rather weak) bounds. 

B. Critical behavior 

The numerical simulations have been performed in net- 
works that are built on a finite periodic one-dimensional 
lattice, except for seed simulations started from the sur- 
face site, where networks built on open chains have been 
used. 

The size (i.e. the number of sites) L of aperiodic net- 
works is not arbitrary but it is given by the possible 
lengths of finite strings of the corresponding aperiodic 
sequence (24|. The system sizes that we have typically 
used in the simulations are shown in Table I and II. The 
simulation time was typically 2 18 - 2 22 MC steps and the 
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ID 


silver- mean 


Fibonacci 


tripling 




1 


0.7864... 


0.7610... 


0.6309... 


drw 


2 


1.7864... 


1.7610... 


1.6309... 


L 




665858 


1346270 


1594324 


ln(6) 




1.7627... 


1.4436... 


1.0986... 


In(r) 




2.5(1) 


2.0(1) 


1.4(1) 


ln(-r) 
ln(6) 




1.4(1) 


1.4(1) 


1.3(1) 


A c 


3.29785 


3.0831(4) 


3.0146(2) 


2.79926(5) 


Vo 


0.31368 


0.261(3) 


0.253(1) 


0.192(4) 


Vi 




0.386(2) 


0.385(2) 


0.386(2) 


»7av 




0.310(3) 


0.307(3) 


0.302(3) 


V« 


0.04998(2) 


-0.080(2) 


-0.067(3) 


-0.016(2) 


So 


0.15947 


0.249(3) 


0.264(3) 


0.367(5) 


S t 




0.124(2) 


0.131(2) 


0.172(2) 


Sav 




0.201(3) 


0.207(3) 


0.253(3) 


5s 


0.42317(2) 


0.590(2) 


0.583(2) 


0.570(2) 


2/z 


1.26523 


1.424(4) 


1.445(3) 


1.618(6) 


2/z t 




1.429(3) 


1.453(3) 


1.625(3) 


2/z s 


1.26523 


1.426(5) 


1.451(5) 


1.624(4) 


a 


0.15947 


0.201(1) 


0.205(1) 


0.250(1) 





k=2 tripling 


k—2 silver-mean 


Hanoi 




0.6309... 


0.6232... 


0.5 


drw 


1.4650... 


1.4575... 


1.3057... 


L 


1594323 


1607521 


4194303 


ln(6) 


1.0986... 


1.7627... 


0.6931... 


ln(r) 


1.4(1) 


2.1(1) 


- 


ln(-r) 
ln(f>) 


1.3(1) 


1.2(1) 


- 


A c 


2.31269(3) 


2.28979(3) 


2.18432(4) 


Vo 


0.117(2) 


0.114(4) 


0.12(1) 


Vi 


0.303(5) 


0.308(2) 


0.33(1) 




0.296(1) 


0.294(2) 


0.28(1) 


Vs 


0.119(2) 


0.117(2) 


0.07(3) 


So 


0.448(4) 


0.453(3) 


0.49(1) 


Si 


0.261(2) 


0.259(2) 


0.28(1) 


Sav 


0.268(2) 


0.271(2) 


0.33(1) 


Ss 


0.446(2) 


0.448(2) 


0.54(1) 


2/zo 


1.665(2) 


1.673(2) 


1.90(1) 


2/zi 


1.665(2) 


1.675(3) 


1.90(1) 


2/z s 


1.667(4) 


1.674(2) 


1.89(1) 


a 


0.269(1) 


0.271(1) 


0.33(1) 



TABLE I: Various properties and estimated critical exponents 
of k = 1 networks. The known ID exponents are from Refs 
1] and 0. 



TABLE II: Various properties and estimated critical expo- 
nents of k = 2 networks. 



averaging has been performed over 10 6 independent runs. 
In case of seed simulations, the applied system sizes were 
large enough compared to the size of growing clusters 
such that the system can be regarded as infinite. 

In order to estimate the critical infection rate A c and 
to keep track corrections to scaling more clearly we have 
monitored the effective exponent a e s(t) defined by the 
local slope 



dln(i) 



(8) 



This kind of analysis helped to estimate the other dy- 
namical exponents as well. However, the pres- 
ence of log-periodic oscillations made the determination 
of the critical point rather difficult, since they distort 
the monotonicity of the functions. Without these modu- 
lations the effective exponents show upward(downward) 
curvature above(below) the transition point, respectively 
(see Fig. S]). In what follows, we shall illustrate the crit- 
ical behavior of the model mainly for the particular case 
of the k = 1 Fibonacci network; the behavior of the pro- 
cess on the other networks is qualitatively similar and the 
corresponding quantitative data can be found in Table I 
and II. 

First we located the critical point by measuring dy- 
namical quantities and calculating the effective expo- 
nents. In case of the k = 1 Fibonacci network, the aver- 
age number of particles originating from the central site 



increases algebraically, superimposed with log-periodic 
oscillations as can be seen in Fig. 0J The critical point 



2.4 
2.2 

2 

CO 

S 1.8 

O 

1.6 
1.4 
1.2 

1 



0.28 




10° 1 o 2 1 1 0" 10' 



1 2 3 4 5 6 7 

ln(t) 



9 10 11 12 13 



FIG. 4: Time-dependence of the particle number No(t) in the 
k = 1 Fibonacci network for A = 3.013, 3.014, 3.0145, 3.016 
(from bottom to top). Log-periodic modulations at criticality 
with a period lnr = 2.0(1) can be read-off. The upper inset 
shows the corresponding local slopes; the lower inset shows 
the unsealed data. 

estimated by the local slopes is at A — 3.0146(2) with 
moderate corrections to scaling. The survival probability 
P(t) and the R 2 (t) also show these periodic modulations 
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with the same period lnr = 2.0(1). This is in agreement 
with the expectations for critical systems with discrete 
spatial scale invariance. Such systems remain self-similar 
when lengths are rescaled by a given (non-arbitrary) scale 
factor b that is characteristic for the system. The time- 
dependent quantities in such models are expected to dis- 
play log-periodic oscillations with a temporal period In r 
that is related to the spatial scale factor through 



In t = z In b, 



(9) 



where z is the dynamical exponent of the model [271 ]. 
The numerical values of the spatial scale factors b of the 
networks under study [23|, [24| and the estimated tempo- 
ral periods are given in Table I and II. The data are in 
satisfactory agreement with Eq. Q. 

As aforementioned, we performed simulations with 
three different positions of the initial active seed. Ac- 
cording to the numerical results (see Table I and II), in 
the three cases the exponents rj m and S m are different, 
whereas z m and the sum r\ m + S m are independent of the 
initial position m [35[ . The latter means that the growth 
rate of the number of active sites averaged in surviving 
samples, i.e. N m (t)/P m (t), is not influenced by the loca- 
tion m of the initial seed. 

As opposed to this, the survival probability is sensi- 
tive to this circumstance. We have found that So > <5;, 
which is intuitively obvious and can be explained as fol- 
lows. The survival probability is greatly influenced by 
the relative position of the growing cluster with respect 
to the long edges as the active region reaches them. If 
the process starts from the central site, the cluster grows 
typically symmetrically around the central site. Since the 
longer edges are also located symmetrically around the 
central site, the growing cluster overlaps with itself as in 
a finite system with periodic boundary conditions. This, 
however, decreases the survival probability because the 
rate of unsuccessful activation attempts is higher in an 
overlapping front. Thus in this case, the long edges of 
the network are not utilized in a favorable way from the 
point of view of the survival of the process. Apparently, 
for initial sites which have an environment identical to 
that of the central one but only within a finite radius, 
the growth is described by the same exponents until the 
cluster is within this radius. 

Contrary to this, when the process is started at a site 
with a long edge, the infection is transferred immediately 
to a far away place and the two clusters spreading out 
from the two sites connected by the long edge do not hin- 
der each other. Of course, in case of a long edge of length 
I, the two advancing fronts meet after time t ~ 0{l z ) and 
the exponents rji and Si describe the spreading dynamics 
only below this time scale. At this time scale, a crossover 
occurs to a region with a faster decaying survival proba- 
bility. In fact, for a typical initial seed location, the dy- 
namical quantities suffer a series of crossovers, depending 
on the relative position of the initial site with respect to 
the longer and longer edges the cluster hits. 

On the grounds of the above argumentations, the most 



favorable initial site for the cluster survival is the site 
at the longest edge, while, disregarding the surface site, 
the most unfavorable site is the central one. Therefore, 
the decay of the average survival probability must be 
bounded by the surviving probabilities in the above two 
extremal situations and we expect Sq > S av > Si to hold. 
According to our numerical results, these relations are 
valid indeed. 

Starting from a fully occupied initial state, the decay 
of the density at the critical point is characterized by 
the exponent a as shown in Figj5] As can be seen from 




FIG. 5: Density decay in the k — 1 Fibonacci network for A = 
3.0145,3.016 (from bottom to top). Log-periodic oscillations 
can be seen as for No(i). The inset shows the unsealed data. 

the estimated exponents given in the tables, in general, 
S m i= ot, in contrast with the CP on regular lattices, 
where S = a holds. As a consequence, the hyper-scaling 
law of DP does not hold on the networks under study, i.e. 
2<5 m + Vm 7^ d/z with d = 1. But the set of critical expo- 
nents are compatible with the generalized hyper-scaling- 
law [H: 



S m + a + r) m = d/z . 



(10) 



with d = 1. We thus conclude that the rapidity reversal 
symmetry, which is necessary to the validity of the hyper- 
scaling law is broken. In other words, P(t) and p{t) scale 
in a different way, which is a consequence of the presence 
of long edges that render the system inhomogeneous in 
space. The inhomogeneity of the system is illustrated 
in Fig. where the local particle densities (ni(t)) are 
plotted at different times against the distance I measured 
from the center. As can be seen, the regularly arranged 
long edges induce modulations in the profiles. As the 
number of infected sites grows at criticality as given in 
Eq. 2] and the active sites are concentrated in a region 
of size ~ t 1 ^ , we expect the density profiles to have 
the following scaling form: 



(11) 
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This has been plotted in Fig. [6] As can be seen, the 
number of peaks in the profiles is increasing with t and, as 
a consequence, the scaling function p(x) is non-smooth. 
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FIG. 6: Upper figure: Average occupation numbers {m(t)) 
of sites at different times, plotted against the distance I mea- 
sured from the center. The logarithms of times are integer 
multiples of the estimated temporal period lnr. The process 
was started from the central sites of the k — 1 tripling net- 
work and the averaging has been performed over I0 7 runs. 
Lower figure: Scaling plot with the same data. 

Nevertheless, the numerical results suggest that S av — 
a is valid within error margin, thus for the average 
cluster-spreading exponents the hyper-scaling of DP is 
fulfilled, i.e. 



of critical exponents, which depend on the underlying 
network, can be found in Table I and II. 



2.2 




FIG. 7: Time-dependence of the particle number No(t) in the 
k — 1 silver-mean network for A = 3.081, 3.082, 3.083 (from 
bottom to top). The period of oscillations is lnr = 2.5(1). 
The inset shows the effective exponents. 




In(t) 



2S av + Tj av = 1/z. (12) 

This indicates that the breaking of rapidity reversal sym- 
metry is indeed related to the presence of the spatial inho- 
mogeneities in the system. This is in agreement with field 
theory of directed percolation with long-range spreading 
[29l ] . where the rapidity reversal symmetry persists. 

Similar to the k — 1 Fibonacci network, we performed 
the above analysis for the k — 1,2 silver-mean and 
tripling networks and for the k — 2 Hanoi-tower net- 
work (see Figs. I7l8l andl5|). The contact process on these 
networks exhibits the same qualitative features as on the 
k = 1 Fibonacci network except that, for the Hanoi-tower 
network, log-periodic oscillations cannot be observed pre- 
sumably owing to their small amplitude. The estimates 



FIG. 8: Time-dependence of the survival probability Po(t) in 
the k = 2 silver-mean network for A = 2.2896, 2.28979, 2.29 
(from bottom to top). The period of oscillations is lnr = 
2.1(1). The inset shows the effective exponents. 



V. DISCUSSION 

First, we give a brief summary of the results obtained 
so far. In each network, a phase transition between ac- 
tive and inactive phases can be identified at some finite 
value of the control parameter by inspecting dynamical 
properties such as the time dependence of the number 
of active sites, the survival probability and the spread. 
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FIG. 9: Time-dependence of the spread Ro(t) in the k — 2 
silver-mean network for A = 2.2896, 2.28979, 2.29 (from bot- 
tom to top). The inset shows the unsealed data. 



At the absorbing phase transition, conventional power- 
law dependence of the above quantities can be observed 
apart from log-periodic oscillations that are related to the 
discrete scale invariance of the underlying networks. In a 
given network, the critical exponents r\ m and 5 m depend 
on the location of the initial seed. The cluster exponents 
satisfy the generalized hyper-scaling relation. However, 
in case of averaging over runs started at different initial 
seed coordinates, the hyper-scaling relation of DP holds, 
too. This means that the rapidity reversal symmetry is 
broken as a consequence of the spatial inhomogeneity. 
The dynamical critical exponents of the phase transition 
differ from that of the one-dimensional DP universality 
class and are found to be characteristic for the underlying 
network. 

We close this work by discussing the relation between 
the critical exponents of the CP and the shortest-path 
dimension (or random-walk dimension) of the underly- 
ing network. The measured critical exponents of the six 
network models are plotted against the random-walk di- 
mension of the corresponding network in Fig. IIP! For 



FIG. 10: The critical exponents r/ av (•), S av = a (M) and z 

(♦) of the six studied systems and the one-dimensional CP 
plotted against d rw of the underlying network. 



decreasing d lw the critical exponents rj av and 8 av = a 
move towards the mean- field values of DP (tjmf = 0, 
Smf = 1) almost monotonically. These exponents lie in 
between the corresponding values of the one-dimensional 
and the two-dimensional DP universality classes. Con- 
trary to this, the dynamical exponent z does not move 
toward the mean field value zmf — 2, but decreases with 
decreasing d rw . Thus, it moves parallel with the dynam- 
ical exponent of the random walk. 

We recall that the k = 1 silver mean and k = 1 Fi- 
bonacci network are the first two members of a family 
of networks, which are defined by the inflation rule in 
Eq. (TT]) , and which are parameterized by an integer n. For 
these networks, it is known that d m i n — ► 1 and d rw — > 2 
when n — ► oo, i.e. in this limit, the characteristics of the 



regular one-dimensional lattice are recovered [24J . There- 
fore, we conjecture that the critical exponents of CP de- 
fined on these networks approach the one-dimensional 
DP values without limits as n — > oo. In the opposite 
limit, i.e. when d rw — > 1 (and d m i n — > 0), it is an open 
question what are the limiting values of the critical ex- 
ponents of the CP. 

Finally, we mention that the above observations open 
up the possibility to design networks on which dynamical 
processes evolve in a prescribed way. The optimization of 
spreading or transport processes in networks is of great 
practical relevance pll l30l l3l| . Besides the tuning of 
degree-distribution or the dependence of transition rates 
on the de gree s of sites in non-regular (such as scale-free) 
networks [32l l33l . l34j , the networks studied in this work 
offer an alternative way of controlling critical dynamics 
by means of choosing appropriate sets of long links and 
the findings obtained here may provide ideas in optimiza- 
tion problems. 
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